Frustrated Polyelectrolyte Bundles and T = Josephson-Junction Arrays 
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We establish a one-to-one mapping between a model for hexagonal polyelectrolyte bundles and a 
model for two-dimensional, frustrated Josephson-junction arrays. We find that the T = insulator- 
to-superconductor transition of the quantum system corresponds to a continuous liquid-to-solid 
transition of the condensed charge in the finite temperature classical system. We find that the role 
of the vector potential in the quantum system is played by elastic strain in the classical system. Ex- 
ploiting this correspondence we show that the transition is accompanied by a spontaneous breaking 
of a discrete symmetry associated with the chiral patterning of the array and that at the transition 
the polyelectrolyte bundle adopts a universal response to shear. 
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Negatively charged DNA molecules in aqueous solu- 
tion condense into dense hexagonal bundles (or tori) in 
the presence of low concentrations of positively charged 

multivalent ions, or "counter-ions" U 13, Ej. This 

property, which has interesting applications in biology 
and which is exhibited as well by other highly charged 
biopolymers, has attracted significant theoretical inter- 
est because it is in clear disagreement with the classical 
Poisson-Boltzmann mean-field theory of aqueous electro- 
statics. Numerical Q and analytical 0, Q studies of 
the so-called "primitive model" -where the biopolymer is 
treated as a linear, charged rod and the counter-ions as 
point charges-indicate that the attraction between two 
adjacent rods is in fact a fundamental feature of aqueous 
electrostatics. The interaction, which has a short range, 
is produced by "out-of-phase" correlations between or- 
dered counter-ion arrays on the two rods. The purpose 
of this letter is to show that the problem of the finite 
temperature counter-ion freezing, or "Wigner crystalliza- 
tion" H , of a hexagonal polyelectrolyte bundle can be 
mapped onto a T = quantum phase transition, and fur- 
ther, that this mapping leads to surprising predictions for 
the properties of polyelectrolyte bundles. 

To demonstrate this claim, we begin with a single rod 
with a uniform fixed charge per unit length, —epo, plus 
a neutralizing distribution of mobile polyvalent ions of 
charge, +eZ, condensed onto the rod. The rod is placed 
in a (monovalent) saline aqueous solution with Debye pa- 
rameter, k. In the absence of thermal fluctuations, the 
polyvalent counter-ions form a one-dimensional "lattice" 
with spacing, d = Z/pq. The modulated part of the 
charge density of the rod can be expressed as a sum over 
the reciprocal lattice vectors of the Wigner crystal: 
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Here, G — 2n/d is the reciprocal lattice vector, \p n \ an 
amplitude determined by the charge distribution of a sin- 
gle counter-ion, and 4>{z) a phase variable restricted to 
the interval [0,27r]. Specifically, (f>(z)/G is the local dis- 



placement of the of the counter-ion lattice with respect 
to the rod. The phase variable obeys an effective Hamil- 
tonian, 



(2) 



where CG 2 is the one-dimensional, T = compression 
modulus of the lattice-roughly proportional to \pi | 2 e ^ Kd - 
and L is the length of the rods. Due to phase fluctuations 
at finite temperatures, the contribution from the n th term 
in eq. to the thermal expectation value (Sp(z)} van- 
ishes as e~ n L /I 3C and the distribution has only liquid- 
like correlations, (5p(0)Sp(z)) ~ e ~\ z \/^ cos Gz, with 




FIG. 1: The two degenerate, "chiral patterned" ground states 
of the antiferromagnetic XY model on a triangular lattice 
are shown in (a). The helical ordering around a triangular 
plaquette of the bundle lattice corresponding to one of these 
ground states is shown in (b). The dark black lines depict 
the polymer backbone while the spheres depict peaks in the 
condensed charge density. 
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£ 2 = 2/3C. Below, we include only the lowest order terms, 
n = ±1. 

We now consider a large hexagonal bundle of rods with 
rod-rod spacing, a, which is maintained by a "bare" ra- 
dial potential that does not depend on the counter-ion 
degrees of freedom discussed here. The bundle is assumed 
to have the elastic properties of a columnar liquid crys- 
tal in the absence of charge ordering. The modulated 
charged density of individual rods given by 0) will then 
be coupled by a Coulomb interaction between the rods, 

U int = J dzidzj 5pi(zi)V(\ri - r 2 \)Spj(zj) 

~ Eij J dz cos (<pi(z) - 02 (z)) ■ ( 3 ) 

The coupling strength, E^, is roughly proportional to 
\pi\ 2 e Kaij with dij the spacing between rods i and j and 
k = V k 2 + G 2 is an effective screening parameter. This 
expression is valid in the limit of slow phase variations, 
where d z (f>i <C G. The effective Hamiltonian of the phase 
fluctuations of the bundle is then, 




where the sum Y]/jj\ is over nearest-neighbor pairs. If 
the bundle is not elastically deformed, then Aij — tt and 
the Coulomb coupling energy is minimized by the afore- 
mentioned out-of-phase condition, <f>i(z) — <fij{z) = n. In 
a hexagonal bundle this condition cannot be satisfied for 
any triplet of adjacent rods (see Fig. QJ; thus, the Wigner 
crystal is frustrated. 

The energy cost of deforming the polyelectrolyte bun- 
dle in the absence of the counter-ion modulation will be 
expressed as a generalized elastic strain energy, 

F[e] = \j d\d?r'e l3 {v)G% kl {v-v')e kl {v') , (5) 

where ey = (diUj + djUi)/2 is the strain tensor and 
u(r) is the displacement field for the bundle. For a 
conventional hexagonal solid, the strain correlation func- 
tion, G?- W (r), would have the form Cfj kl S(r), with C°- fci 
the tensor elastic constants having five independent en- 
tries 0- If the bare rods are uniform-as we assume- 
then entries such as C xzxz = Cy Zyz that describe resis- 
tance to relative sliding of the rods along the I direc- 
tion must vanish. That case corresponds to the colum- 
nar liquid crystal where the elastic energy contains a 
term, (1/2) J d 3 rK (d 2 u x /dz 2S ) 2 , where K is the bending 
stiffness. The Fourier transform of the strain correlation 
function is therefore G xzxz (c[±, q z ) ~ Kq 2 . The develop- 
ment of the counter-ion modulation can be viewed as a 
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FIG. 2: Two volume-preserving, shear deformations which 
lead to the relative sliding of the neighboring rods (equivalent 
up to a pure rotation). 

phase transition between a columnar liquid crystal and 
a hexagonal crystal. Because the phase variable, </>i(z), 
describes sliding of the counter-ion lattice relative to the 
jth roc ^ -j-jjg optimal phase difference between rods i and 
j will be altered by any shear strain that produces either 
a relative sliding or a simultaneous rotation (see Fig. |2J) 
of the rods, 

A i:j (z) = 7T + 2G / dr ±k e kz (r ± ,z) . (6) 

J i 

The integral here is to be taken along a path in the x-y 
plane that connects the two rods. 

For the case when Ay = tt and only nearest-neighbor 
rod couplings are considered, the partition function, Z — 
f[d(/)i]e~P F W i \ can be mapped onto the partition func- 
tion of a T = quantum, many-body system, namely a 
fully frustrated two-dimensional Josephson-junction ar- 
ray with a / = 1/2 of a flux quantum per elementary pla- 
quette, and with Coulomb interaction restriced to capac- 
itive charging of the individual grains [l(3 | . The partition 
function of such an array can be expressed as the path in- 
tegral of exp { — S[(j>i(T)]/h} over all phase configurations 
<fri i T ) with r an imaginary-time variable and the classical 
action, h~ 1 S[4>i(T)] equal to j3F[<j>i(z)] with z — > r and 
L — > hp q , where (3 q corresponds to the temperature of the 
quantum analog system. Here, we see precisely that the 
infinite rod limit, L — > oo, of the classical bundle system 
maps onto the T — > limit of the quantum superconduct- 
ing array. In the quantum system 4>i(r) represents the 
phase of the superconducting wavefuction conjugate to 
the number operator (rij = —id/d(f)i) for Coop er-p aired 
electrons on the superconducting grain at site i |ll| . The 
phase stiffness, /3C, in eq. corresponds to the in- 
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verse of the capacitive charging energy of a grain, fiE' 1 , 
and the nearest-neighbor Coulomb coupling, /?.Ey, cor- 
responds to the Josephson inter-grain coupling energy, 
h~ 1 Ej. Table [I] gives a full list of corresponding quanti- 
ties. 

As a function of the ratio of charging and Josephson 
energies, a = E c /Ej, such an array undergoes, with de- 
creasing a, a continuous quantum phase transition at 
a c = 6 from an insulating to a superconducting phase 
which has a non-zero expectation value for the operator, 
e z0i ^phe orc i er parameter, ^ -related to (e Ic ^)-for 

this transition vanishes near the critical point as \a— a c \° 
and the superfluid density, p si as \a — a c \ v with critical 
exponents j3 — 0.25(2) and v = 0.50(2) computed for the 
d = 3, n = 2 stacked triangular XY anti-ferromagnet 
[l2|. These results have interesting implications in the 
context of polyelectrolyte bundles. The control param- 
eter, a, corresponds to (f3 2 CE i j)~ 1 and (e 4< ^) to the 
expectation value of the modulated charge density of 
the Wigner crystal, (Spi(z)). This mapping evidently 
predicts a continuous freezing transition of the counter- 
ions as a function of temperature, even though three- 
dimensional freezing transitions are in general first-order. 
The critical exponents are, in fact, close to those of a tri- 
critical point [l2j , which does indicate the proximity of 
a first-order transition. 

A further surprise is due to the fact that an / = 1/2 
Josephson-junction array in the superconducting phase 
also has broken discrete symmetry with an associated 
Z2 x £7(1) order parameter [Jjj. This is due, as shown 
in Fig. ^ to the fact that there are two degenerate sets 
of chiral-patterned ground states for a d = 2 triangular 
XY anti-ferromagnet, with the phase winding around 
a particular plaquette in opposite senses and each spin 
27r/3 out of phase with its neighbors 0|. In particular, 
the order parameter, 'J, can be constructed from a pair of 
complex numbers, ip + and ip~ , related to the expectation 
value of e l ^ i by, 

(e 4 ^) cx V + e +iq *- Xl + V^e- 4q *- x * , (7) 

where q* = (47r/3a)cc, is the wavevector associated with 



TABLE I: Correspondence between the d — 3, classical poly- 
electrolyte system (left) and the d = 2+1, quantum frustrated 
Josephson array system (right). 



free-energy functional, /3F[<j>] 
height along polymer, z 
longitudinal phonon 

(or phase) stiffness, /3C 
electrostatic, inter-rod 

coupling, pEij 
shear strain, 2Ge± z 
elastic response, G xzxz (q z ) 



imaginary-time action, h~ 1 S[ 
imaginary time, t = it 
inverse grain charging 

energy, hE^ 1 
Josephson inter-grain 

coupling, h~ Ej 
vector potential, (27r/$o)a± 
current response, — C xx {— icu) 



the helical winding of the phase. For a < a c the dis- 
crete symmetry is broken and the order parameter spon- 
taneously aligns along the + or — direction. This indi- 
cates that in the charged-ordered phase the modulated 
density, (<5<fo), adopts a locally helical configuration (see 
Fig. ^) , where the winding on neighboring plaquettes has 
the opposite handedness. 

The electromagnetic properties of a Josephson array, 
such as conductivity and diamagnetic susceptibility, are 
determined by the response of the array to an applied 
vector potential. Specifically, the conductivity can be 
expressed by the Kubo relation as <t{lu) = lu\C xx {(\ = 
0, u)/u) where 



G Lj (q^,-iuj) 



5 2 \nZ 



5oi(r, T)8aj(0, 0) 



, (8) 



is the analytic continuation of the current-current cor- 
relation function to imaginary frequencies |14| , a(r^ , t) 
is an infinitesimal vector gauge field added to an exter- 
nally applied vector potential, i.e. A — > A + a. In the 
context of polyelectrolyte bundles, the externally applied 
vector potential would correspond to the = n con- 
dition for the undeformed bundles while, according to 
eq. JiU , the effect of applying an infinitesimal gauge field 
(27r/<f>o)afc(r^, t) on Ay corresponds to that of intro- 
ducing a shear strain, (2G)e z fc(r^, z). It follows from the 
definition of the elastic strain energy and eq. © that the 
finite-temperature strain correlation function equals 



G ijk i (r) 



Se lj (r)Se k i(0) 



(9) 



Comparison of eqs. JBJl and shows that the 
Fourier transform G xzxz (q±, q z ) corresponds to minus 
the C xx (q± 1 —iuj) component of the current-current cor- 
relation function. Thus, we can make use of established 
results for the current-current correlation function of the 
Josephson array to deduce certain elastic properties of a 
polyelectrolyte bundle. 

In the insulating phase (a > a c ) of the Josephson ar- 
ray, the low-frequency conductivity is linearly propor- 
tional to the frequency and C xx (—ioj) cx — w 2 £(a), where 
is the correlation length which diverges at the crit- 
ical point as \a — a c \~ v ' ■ For the polyelectrolyte bun- 
dle, this means that in the molten phase with no re- 
sponse to uniform sliding deformations (or C xzxz = 0), 
G xzxz (q±,q z ) cx (k B TG 2 )q 2 £,(a). Hence, we have a 
columnar liquid crystal with a bending stiffness, K, of 
the bundle which diverges at the critical point as the 
correlation length, £(a). Turning next to the supercon- 
ducting phase (a < a c ) of the Josephson array, since the 
zero-frequency conductivity is infinite the zero-frequency 
current-current correlation function is proportional to the 
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superfluid density, so that C xx (c[± — 0,co = 0) oc £ 1 (a) 
1 1 51 ] - Therefore, the shear modulus of the polyelectrolyte 
bundle, C xzxz = G x zxz (l = 0) cx (fc B TG 2 )r is fi- 
nite in the low-temperature, charge-ordered phase and 
vanishes at the critical point. This is the modulus corre- 
sponding to the uniform elastic deformations pictured in 
Fig. [21 Finally, right at the critical point (a = a r ) 
the Josephson array is predicted to be metallic 
i.e. to have a finite zero-frequency conductivity with 
C XI (qi, — iu) — <J*\tu\. Moreover, the critical conduc- 
tance is expected to be a universal quantity equal to, 
The corresponding strain correlation func- 
tion, G xzxz (q±,q z ) oc (k B TG 2 )\q z \, would belong to an 
unusual mesophase with an infinite bending stiffness, 
unlike conventional columnar liquid crystals, yet unlike 
hexagonal crystals which have infinite bending stiffness, 
it has no broken translational symmetry along the z di- 
rection and hence no resistance to sliding deformations. 

Another interesting correspondence involves the topo- 
logical excitations associated with the gauge field. In 
the Josephson-j unction array a closed-line integral of the 
gauge field 

Q-) /dr x .a = 0,±l,±2,... (10) 

corresponds to one or more flux quanta added or removed 
from the system. An added flux quantum is a stable topo- 
logical excitation of the superconducting phase but not of 
the normal phase. For a polyelectrolyte bundle the anal- 
ogous quantity, the closed line integral § drk(du z /drk), 
corresponds the the Burgers vector of a screw dislocation 
long the z direction Screw dislocations are indeed 
well-defined topological excitations of hexagonal crystals 
but not of columnar liquid crystals. The energetic cost 
per unit length of the screw dislocation is proportional 
to the C xzxz = C yzyz shear elastic constants, which van- 
ish at the critical point as £(a) _1 . This suggests that 
one could view the counter-ion freezing transition as an 
unbinding transition for screw dislocations loops. 

In more realistic models of biopolymers important fea- 
tures such as the excluded core region, the chirality of the 
rods, and non-pairwise rod-rod interactions would have 
to be included. A preliminary study indicates that mod- 
els incorporating such features, while more complex, still 
display a counter-ion freezing transition closely similar 



to the one discussed in this letter. An important general 
conclusion of this work we expect to hold in more realis- 
tic models is that the elastic constants of polyelectrolyte 
bundles are of fundamental interest for understanding the 
phase behavior of the counter-ions. Interestingly, a twist 
elastic deformation due to counter-ion condensation has 
been observed for the case of F-Actin bundles 0|. A 
study of the freezing transition for more complex, chiral 
geometries appropriation for such cases is in progress. 
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